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ABSTRACT 



We investigate the shapes of 7— ray pulsar light curves using 3D pulsar mag- 
netosphere models of finite conductivity. These models, covering the entire spec- 
trum of solutions between vacuum and force-free magnetospheres, for the first 
time afford mapping the GeV emission of more realistic, dissipative pulsar mag- 
netospheres. To this end we generate model light curves following two different 
approaches: (a) We employ the emission patterns of the slot and outer gap models 
in the field geometries of magnetospheres with different conductivity a. (b) We 
define realistic trajectories of radiating particles in magnetospheres of different o 
and compute their Lorentz factor under the influence of magnetospheric electric 
fields and curvature radiation-reaction; with these at hand we then calculate the 
emitted radiation intensity. The light curves resulting from these prescriptions 
are quite sensitive to the value of a, especially in the second approach. While still 
not self-consistent, these results are a step forward in understanding the physics 
of pulsar 7— radiation. 

Subject headings: pulsars: general — stars: neutron — Gamma rays: stars 



1. INTRODUCTION 



The Fermi Gamma-Ray Space Telescope has had a major impact on our understanding 
of pulsar physics with the discovery of over 100 7— ray pulsars comprising three populations: 
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youn g radio-loud pulsars, young radio-quiet pulsars and millisecond pulsars flAbdo et al. 
2010l ). Studying 7— ray pulsars with such a broad range of underlying physical parameters 
offers an opportunity of deeper understanding the physics underlying the pulsar 7— ray emis- 
sion and magnetic field geometry. A major issue resolved early in the Fermi mission was the 
site of the pulsar high-energy (GeV) emission. The cutoff of the Vela pulsar phas e-averaged 
spectrum measured by the Fermi Large Area Telescope (LAT) (lAbdo et al.l 120091 ) ruled out 
at high significance the super-exponential shape of magnetic pair production attenuation in 
polar cap cascades and established the location of the 7— ray emission and particle acceler- 
ation in the outer magnetosphere. 

Present mo dels for pulsar high-energy emission assume a vacuum retarded dipole (VRD) 
( IDeutschl Il955l ) field geometry, which is expedient but fundamentally inconsistent. Such 
models neverthele ss have had some success in modeling Fermi pulsar light curv es (LCs). 



Outer Gap (OG) (ICheng et al. 



1986 



Romani fc Yadigaroglulll995t iHirotani 1120081 ) and Slot 



Gap (SG) flMuslimov fc HardineJbool models both derive regions of E \ bordering the last 



open field-lines extending to the light-cylinder. The particles accelerating in these 'gaps' 
produce curvature radiation (CR) and inverse- Compton emission, with their Lorentz factors 
limited by CR reaction forces that balance En. The pattern of emission on the sky shows 
caustics that form on the trailing edge of the open field region of each magnetic pole as phase 
shifts due to dipole geometry, time-of-flight and aberration nearly cancel, and p hotons from 
a large range of altitudes arrive in phase (jMorinil Il983l ; iDyks fc Rudakl 120031 1 . Observers 
viewing at angles crossing the caustics will see one or two narrow peaks that resemble the 
7— ray LCs seen by Fermi. Because the GeV emission in these models takes place in a region 
near the pulsar last open field-lines, the shape of the model LCs is sensitive to, and thus a 
good diagnostic of, the geometry of the pulsar magnetosphere near the light-cylinder. 

Fortunately, this can now be addressed in detail, thanks to recent advances in numerical 
simulation of pulsar magnetospheres that model the high-altitude field structure critical to 
the high-energy emission. The global structure of realistic pulsar magnetospheres remains an 



ometries of VRD and force-free-electrodynamics (FFE) 


(Contopoulos et al. 


1999 


Spitkovskv 


2006 




Timokhin 


2006; 


Kalapotharakos & Contopoulos 


2009 


). The effects of acceleration 


fields 


Hirotani 


2007. 


20081) and open-zone currents (] 


Romani & Watters 


20HT 


on the LC 



have been explored, but these models are not fully self-consistent. In all these models, the 
sweepback of the magnetic field-lines near the light-cylinder, due to retardation and currents, 
causes an offset of the polar cap (PC) (and of the entir e magnetosphere) i n the direction 
opposite to the rotation, which can affect the 7— ray LCs. iBai fc Spitkovskyl ( 120101 ) modeled 
7— ray LCs in F FE field geometry injecting pho t ons al ong tangents to the field direction of 
the separatrix. IContopoulos fe Kalapotharakos! ( 120101 ) injected photons only in regions of 
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the FFE magnetosphere with J/pc = 1, where J is the current and p is the local charge 
density, as in these regions the electron velociti es are expected to be sufficiently close to c to 
lead to GeV photon production. More recently, Harding et al. (j2011 ) assumed SG geometry 



to produce model LCs for the VRD and FFE magnetospheres, concluding that the VRD 
geometry provides better fits to the observed LCs than the FFE geometry. The FFE mod- 
els present larger field-line sweepback and consequently the corresponding LCs have larger 
phase-lags relative to the radio pulse which is not consistent with that of observed Fermi 
LCs. 

Most recently, resistive magnetosphere models have appeared in the literature (Kalapotharakos 



et al. 2012 (K12, hereafter) ; iLi et al.ll2012l ). that drop the ideal-MHD requirement in favor 
of an Ohms' Law that relates the current to the E and B fields through a finite conductiv- 
ity. These simulations reveal a range of magnetic field structures, current distributions and 
spin-down power that lie between the VRD and the FFE solution. Most importantly, mod- 
els of finite conductivity possess regions of E\\ which are potential locations of high-energy 
emission. 

In this Letter, we explore for the first time the high-energy emission that is generated 
in resistive magnetospheres, using their magnetic field structure and E\\ to produce 7— ray 
LCs following two different approaches. First, we assume the emission geometry of the 
SG and OG models and compute the resulting LCs for different values of conductivity to 
compare with those of the VRD and FFE geometries. In the second approach, we define 
approximate particle trajectories, and we calculate the corresponding energies (including 
radiation losses) and CR emission. While neither of these approaches is self-consistent, in 
that particle motions should affect the fields, they are an important first step in relating the 
field structure and acceleration dictated by global magnetosphere solutions to observations. 



2. LIGHT CURVE MODELING 

2.1. Pulsar Magnetosphere Models with Finite Conductivity 

Dissipative magnetospheres are necessary for modeling pulsar LCs considering that nei- 
ther the VRD nor the FFE solutions are compatible with the emission of radiation: The VRD 
solutions provide maximum accelerating field, En, but no charges (p = 0), while the FFE 
solutions have a sufficiently large number of particles to guarantee the nulling of En . Herein, 
we use the dissipative solutions presented recently by K12, which use a phenomenological 
conductivity a {in lieu of microphysical processes) to relate the current density J to the 
fields E, B. Although these solutions are still not self-consistent, they are an improvement 
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over those of VRD and FFE because, besides the global field geometry, they also provide the 
distribution of Eh, a quantity necessary to compute the acceleration of radiating charges. 

In the next sections we present model LCs using the structure of magnetic and electric 
fields provided by dissipative solutions of the perpendicular rotator (a = 90°) which span 
the entire solution space from VRD to FFE. These solutions are presented in detail in K12 
and have been produced adopting a very simple prescription for the current density 

E x B 

J = cp— ^- + aEy (1) 

In this case the current density consists of two components, namely a drift current and a 
component parallel to the magnetic field. While some simulations explored models with a 
spatially dependent a, in this paper we will use only those with constant a. As a goes from 
to oo the corresponding solution ranges from VRD to FFE. The solutions we consider here 
correspond to er ~ 0.08, 1.5 and 240, where O is the angular frequency of the star. The field 
structure of the solutions for a ~ 1.50 and o th 240 are shown in the last rows of figures 4 
and 3 of K12, respectively. 



2.2. Geometric Approach 



A simple method of generating pulsar LC s, used in many previous studies (e.g lDyks fc Harding 



20041 ; IWatters et al.ll2009t Venter et al.ll2009l ). adopts the geometry of physical emission mod- 
els that have computed the shape of accelerator gaps with assumed field structure and sources 
of charge. This method can directly compare LCs from magnetospheres having finite con- 
ductivity with LCs in VRD and FFE magnetospheres. To explore how the magnetic field 
structure and offset PCs influence 7— ray pulsar LCs, we have generated model LCs using 



a geometrical version of the SG and OG models (e.g. iDyks fc Harding! l2004t IWatters et al. 
2009c Venter et al.l 120091 ) . The SG has its origin in PC pair cascades that screen the accel- 
erating parallel electr ic field . En o ver most of the open field except in narrow gaps along the 
last open field-lines (lAronsl 1 1 98 31 1 . The electrons accelerate and radiate from the neutron 
star (NS) surface to high-a l titude, and emission occurs throughout the volume of the gap 
( jMuslimov fc Harding! I2004J ; lHarding et al.l 120081 ) . The OG is a vacuum gap that also forms 
adjacent to the last open fi eld- line, above the null charge surface where the corotation charge 



( iGoldreich fc Jurianlll969[ ) changes sign. The gap width is determined by the scr eening of E \ 



by pa ir cascades and emission occurs in a thin region along the gap inner edge (IWang et al. 



2010|). 



Com ponents of the magnet ic field are determined from analytic expressions for the 



VRD (see lDyks fc Harding! 120041 ) and interpolated from numerical simulations for FFE and 
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resistive magnetospheres. The open field boundary on the NS surface (the PC rim) was 
determined via bisection in magnetic colatitude at fixed azimuth values. Open Volume 
radial an d azimuthal Coordinat es (OVC) were then defined inside the open volume of each 



solution ( iDyks fa Harding! l2004j ). We assume that particles travel from the NS surface along 
open field-lines in OVCs and emit radiation tangent to field-lines, uniformly, in the corotating 
frame (CF). We assume that emission is also uniform across a SG of width w = 0.05, as a 
fraction of the open volume, on field-lines originating between r min = 0.95 and r max = 1.0 
on the PC (in units of PC radius) and in a thin layer at r min ~ r max = 0.95 in the OG. 
The minimum and maximum spherical radii of emission are assumed to be the NS surface 
and -R max = 1-2-Rlc; limited by a maximum cylindrical radius of R^l x = 0.95-Rlc, f° r the 
SG, and the null surface and -R max = 1.5-Rlc, limited by R^l x = 0.97-Rlc, for the OG. The 
photon direction is assumed to be tangent to the magnetic field in t he CF, obtained throug h 



a Lorentz transformation from the inertial observer's frame (IOF) (IBai fa Spitkovskyll2010l ). 
The emission direction is then transformed to the IOF (aberration), time-delays are added 
and the emission is accumulated in sky-maps in viewing angle ( and phase with respect to 
the pulsar rotation axis. LCs are then obtained as slices through these maps at constant (. 



2.3. Particle Trajectory Approach 

An altogether different approach to produce LCs that takes approximately into account 
the electric fields present in the specific magnetospheric solution is the following: Since 
we anticipate the velocity of any particle of the magnetosphere to be very close to c, we 
decompose its motion into a drift component and one parallel to the magnetic field. Thus, 
one can write for the particle velocity 

E x B „ B 

u=^c + /c- (2) 

The sign and the absolute value of the factor / is chosen so that the motion of the particle 
be outward and the total modulus of the velocity u be c. This trajectory determ i natio n 



is similar to those of IContopoulos fa Kalapotharakosl ( 120101 ) and IBai fa Spitkovskyl ( l2010l ). 



Under this assumption we calculate the trajectories passing through each magnetospheric 
point inside a central cube of edge 3i?LC considering that in the dissipative solutions the 
particles do not follow the field-lines in the CF. So, instead of open and closed field-lines 
we determine open and closed trajectories depending on whether they reach (or not) 2i?Lc, 
assuming that the radiating particles follow only open trajectorie^]. This particle trajectory 



1 Our field structures have Ei\ in the closed zone too; we believe that this is due to approximations in our 
field computations and therefore we restrict radiation by particles accelerated only by the 2S|| of the open 



-6 - 



determination allows also the calculation of the local radius of curvature R cr at each point of 
the magnetosphere. Moreover, assuming that each particle starts at the stellar surface with 
a small 7-value (7 < 100) we can calculate its Lorentz factor 7 along its trajectory from 

^7 _ j (lecE\\ 2 g e V 
dt m e c 2 3 R 2 cr m e c 

where q e and m e are the electron charge and rest-mass respectively. This approach allows us 
to have all the information needed to calculate the CR intensity contributed by each point 
of the magnetosphere and so the corresponding sky-maps and the LCs (see Section 2.2). 



3. RESULTS 

Figure 1 shows geometric LCs for SG (left-hand column) and OG (right-hand column) 
emission in VRD (black) and FFE (purple) solutions, and in resistive magnetospheres with 
a = 0.08f2 (red), and 1.5fi (green) for a range of observer angles (. The resistive solutions of 
lowest a are closest to the VRD and indeed, the computed LCs look very similar. However, 
there is a shift to larger phase and a slight broadening of the peaks of the a = 0.08Q LCs. 
As a increases, in both the SG and OG cases, the peaks are shifted even more to larger 
phase and the broadening is more pronounced, with the highest a-value being closest to the 
FFE solutions. Indeed, their LCs look very similar, with only a slight shift in phase of the 
peaks but no further broadening. Overall, there is a distinct progression in the LC shapes as 
conductivity increases. The VRD LCs have the narrowest peaks and the smallest phase-lag 
from the magnetic pole (phase= 0), with peak- width and phase-lag systematically increasing 
with conductivity. 

The LC changes with conductivity result from changes in magnetic field structure. Mag- 
netospheres with low a are "stiffer" and thus have less sweptba ck field-lines and smalle r open 



field volume, while those with higher a have more sweepback (IBai fc Spitkovskyll2010l ). The 
increase of sweepback with a produces a larger shift of the PC which causes the larger phase- 
lag of the LC peaks. The increase in open volume of magnetospheres with large a also causes 
the increase in peak-width, since the gap widths are assumed to be a fraction of the open 
volume. Since peak-width and phase-lag are measurable characteristics of observed 7— ray 
pulsar LCs, this study shows that they could potentially be an important diagnostic of mag- 
netospheric conductivity. Comparison of geometric LCs in VRD and FFE magnetospheres 



has a lready indicated that VRD provides a better match to observed LCs ([Harding et al. 



201 A 



zone. 
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In Figure [2] we plot the LCs for solutions corresponding to three different values of 
a taking into account the physical properties provided by each solution as we describe in 
Section 2.3. The red, green and blue color correspond to a = 0.080, 1.5Q, 240, respectively. 
For these LCs we assume that emission occurs only along all open trajectories (i.e. those 
that reach at least up to 2i?Lc)- The emission is considered to be due to CR and is always 
proportional to 7 4 i?~ 2 . The 7- value is derived by Eq. [3j The total emissivity can also be 
weighted by the local charge density p (left-hand column of Figure [2]). The general feature 
is that the broadest (narrowest) pulses seem to be those corresponding to the middle (high) 
a- value a = 1.50 (a = 240). However, the LCs for a = 240 and for the low (- values exhibit 
high off-pulse emission. This effect decreases when the charge density p weighting is included. 
Moreover, near £ = 90°, the middle a- value (a = 1.50) pulses are weak0, double and narrow. 
The general trend for the o = 0.08O and a = 1.50 solutions is that the phase-lag of the 
pulses with respect to the magnetic poles (phase = 0) decreases with ( although there are 
counter examples. These phase-lags start from values higher than 0.25 (for ( = 45°) and only 
for ( near 90° can reach close to 0.1-0.15. For the high a- value [a = 240) the corresponding 
phase-lags seem to be near the value 0.25 for most of the cases implying a non-monotonic 
behavior with a. 

We checked also the assumption that the emitting particles are mostly those that follow 
the high-voltage trajectoriefl This approach is similar to the geometric one (Section 2.2) 
in the sense that it is supposed that only a part of the magnetosphere contributes to the 
emission. The difference here is that the active region is traced by the high-voltage trajec- 
tories. Figure |3] is similar to Figure [2] but only the 10% of the highest-voltage trajectories 
are considered to emit. In this case we observe narrower pulses. However, only for the 
middle a- value (or = 1.50) we have pulses corresponding to small phase-lags. We note also 
that the low a-value (a = 0.08O) is observable only from a relatively narrow range of ( 
(C ^45° -60°). 

We derived also the LCs considering each point's emissivity oc p s 7 4 i?~ 2 where p s is the 
charge density on the star surface at the point where each particle starts its journey. In this 
case the results are similar to those of the right-hand columns of Figures [21 [3j 

In Figure H] we plot the points in the 3D magnetosphere that produce the pulses shown 
in the right-hand column of Figure [2l In each of these cases we have identified the phases 
of the observed pulses for all ( > 30° and we located the points that contribute over 90% 
of the corresponding emission. The emissivity of each point is represented by the indicated 



2 This is not shown in the normalized LCs of Figure [5] but only in the corresponding sky-maps. 
3 The voltage is given by J E ■ dl along a trajectory 
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color scale. In these figures we have also plotted the portion of the magnetic field-lines (in 
gray) that contribute to the emission in the SG model presented in Figure [H We see that 
there is always a blob of points over the polar caps contributing to the observed pulses. 
However, as we go toward high a-values the volume of these blobs decreases and new points 
from the outer magnetosphere are added. For a = 0.08Q and a = 1.5Q the inner parts of 
the blobs coincide with only a subset of the SG lines. This subset increases with a and at 
a = 24Q the largest part of the SG lines seem to follow the colored points. Beyond the light- 
cylinder, this region coincides with the current-sheet (Figure 4). This supports the result of 



( IContopoulos &: Kalapotharakosll2010[ ) where a significant part of the emission is produced 
near the current-sheet. The changes of the form of the effective emitting regions, described 
above, make the sky-maps evolution with a complex. This makes the LCs (Figures [2j [3]) 
evolve in a non-monotonic manner with increasing a. However, inspection of the full sky- 
maps shows smooth evolution with increasing conductivity. 

Figures [TH show that the geometric LCs are less sensitive to a at high a- values than 
those of the trajectory approach. This indicates that the En distribution changes more 
significantly than the magnetic field structure at high a. 



4. DISCUSSION 

We have presented a first implementation of dissipative magnetospheres to model pulsar 
emission for direct comparison with the fast accumulating pulsar phenomenology. We have 
concentrated here on the pulsar 7— ray LCs, however our models could also address spectral 
features. With LC modeling alone, we have just scratched the surface of the problem, since 
we have used only one of the prescriptions discussed in K12 and only one inclination angle 
a = 90°. A more complex behavior is expected for a variety of prescriptions and a's. 

The main goal of this paper is to show that more realistic pulsar magnetosphere models 
provide flexibility that allows meaningful constraints on their parameters in direct compari- 
son with observations. While neither of the approaches employed are self-consistent in that 
they take into account the effects of the radiating particles on the magnetosphere itself, we 
find a progression in the LC shape, with peak phase and width increasing with a in the 
SG and OG models. Observed Fermi LAT LCs show an inverse correlation between the 
peak separation A (in LCs with two peaks) a nd phase- lag S of the first peak relative to the 



radio-peak (thought to be near phase = 0) (lAbdo et al.l l2010i ) . The first-peak phase-lags 



(~ 0.17 — 0.2) of high a model LCs are too large to account for the observed S of many 
LAT pulsars with A = 0.5. The LCs computed in the particle trajectory approach can also 
produce narrow pulses, despite the fact that they assume particle emission at every point of 
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the magnetosphere; however, their consistency with observations may require specific values 
of (, a and a. Independently of whether this is true or not, the final calculation of the 
LCs depends also on the modulation by the local number of the emitting particles, which is 
something that can be derived only in fully self-consistent solutions. We plan to move in this 
direction by introducing pair cascade mechanisms and the calculation of the exact particle 
orbits (using the full equations of motion), taking into account both CR and synchrotron 
losses, as well as inverse- Compton radiation. 

We expect that with the dissipative models at hand, comparison of their model LC and 
spectral products with observations will allow us to find specific magnetospheric prescriptions 
and model parameters that provide the working physics of pulsar magnetospheres. All these 
will be the subject of future work. 

AKH acknowledges support from the NASA Astrophysics Theory and Fundamental 
Physics Program and the Fermi Guest Investigator Program. We thank also the referee for 
his/her constructive comments. 
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Slot gap Outer gap 




Phase 



Fig. 1. — Geometric LCs for slot gap (left) and outer gap (right) emission in VRD (black), 
FFE (purple) solutions, and resistive magnetospheres with a = 0.08f2 (red) and cr = 1.5Q 
(green) for pulsar inclination angle a = 90°. 
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phase 

Fig. 2. — The LCs corresponding to the trajectory approach. The assumed emissivity at 
each point of the magnetosphere is indicated in the figure. The red, green and blue colored 
lines correspond to a = 0.08Q, 1.5Q and 24S1, respectively. 
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Fig. 3. — Similar to Figure [2] but only the 10% highest-voltage trajectories are assumed to 
emit. 
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Fig. 4. — The regions of the magnetospheres that produce the peaks of the pulses shown in 
the right-hand column of Figure El The color scale indicates the corresponding emissivity. 
For o = 2AQ a significant part of the emission comes from a region near the current-sheet 
outside the light- cylinder (gray surface in the right-hand panel). 
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